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ABSTRACT 

Carbon  films  of  thickness  down  to  2 nanometers  are  necessary  to  achieve  a storage 
density  of  100  Gbit/in2  in  magnetic  hard  disks.  Reliable  methods  to  measure  the  properties  of 
these  ultrathin  films  still  have  to  be  developed.  We  show  for  the  first  time  that  combining 
Surface  Brillouin  Scattering  (SBS)  and  X-ray  reflectivity  measurements  the  elastic  constants  of 
such  films  can  be  obtained.  Tetrahedral  amorphous  carbon  films  were  deposited  on  silicon,  by  an 
S bend  filtered  cathodic  vacuum  arc,  which  provides  a continuous  coverage  on  large  areas  free  of 
macroparticles.  Films  of  thickness  down  to  2 nm  and  density  of  ~3  g/cm3  were  produced  and 
characterized.  The  dispersion  relations  of  surface  acoustic  waves  are  measured  by  SBS  for  films 
of  different  thickness  and  for  the  bare  substrate.  Waves  can  be  described  by  a continuum  elastic 
model.  Fitting  of  the  dispersion  relations,  computed  for  given  film  properties,  to  the  measured 
dispersion  relations  allows  the  derivation  of  the  elastic  constants.  For  a 8 nm  thick  film  we  find  a 
Young’s  modulus  E around  400  GPa,  with  a shear  modulus  G lying  in  the  130  -210  GPa 
interval.  For  a 4.5  nm  thick  film,  E is  around  240  GPa,  with  G lying  in  the  70-  130  GPa  interval. 
Results  for  even  thinner  films  become  highly  sensitive  to  the  precision  of  the  substrate 
properties,  and  indicate  that  the  above  values  are  lower  bounds.  We  thus  show  that  we  can  grow 
and  characterize  nanometer  size  tetrahedral  amorphous  carbon  films,  which  maintain  their 
density  and  mechanical  properties  down  to  the  nm  range. 

INTRODUCTION 

The  storage  density  of  magnetic  hard  disks  is  steadily  increasing  [1,2],  Technological 
developments  are  required  to  maintain  the  increasing  trend.  Namely,  the  distance  between  the 
recording  head  and  the  magnetic  material  has  to  be  decreased,  and  this  requires  a reduction  of  the 
thickness  of  the  protective  coatings  of  the  disk  and  the  head.  Carbon  coatings  are  presently  the 
standard  choice:  they  provide  protection  against  corrosion,  grafting  sites  for  the  lubricant 
molecules  and  mechanical  protection  [1,2].  The  continuity  of  the  film  is  a crucial  requisite.  The 
achievement  of  a storage  density  of  100  Gbit/in2  requires  a reduction  of  coating  thickness  down 
to  2-5  nanometers.  The  characterization  of  such  thin  films  is  still  an  open  question.  The  elastic 
constants  provide  useful  information:  beside  the  direct  characterization  of  the  mechanical 
behavior,  they  give  indirect  indications  on  the  quality  and  the  continuity  of  the  film.  Although  a 
film  of  few  nanometers  on  a softer  substrate  cannot  provide  a strong  mechanical  protection  [2], 
measurement  of  its  elastic  moduli  provides  however  significant  indications  on  film  quality, 
including  the  fraction  of  sp3  bonding. 

The  stiffness  of  an  amorphous,  and  therefore  mechanically  isotropic,  film  is  fully 
characterized  by  two  independent  parameters,  which  are  typically  taken  among  Young’s 
modulus  E,  shear  modulus  G,  bulk  modulus  B,  Poisson’s  ratio  v or  among  the  elements  Cjj  of  the 
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elastic  tensor.  Any  two  of  them  determine  all  the  others.  Their  measurement  for  ultra-thin  films 
is  arduous.  Indentation  directly  supplies  the  hardness,  and  is  commonly  adopted  to  measure  also 
the  elastic  constants  [3-5].  However  for  film  thicknesses  below  ~1  pm  the  requisite  of 
indentation  depth  significantly  smaller  than  film  thickness  is  met  with  difficulties,  and,  especially 
for  hard  films  on  softer  substrates,  the  behaviour  of  the  substrate  can  affect  the  results.  For 
thicknesses  of  a few  nanometers  indentation  is  no  longer  viable.  The  elastic  properties  must  be 
probed  by  other  methods,  namely  exploiting  acoustic  waves. 

Brillouin  scattering,  i.e.  the  scattering  of  an  electromagnetic  wave  (a  photon)  by  an 
acoustic  wave  (a  long  wavelength  acoustic  phonon)  provides  a convenient,  contact  less  way  to 
measure  the  spectrum  of  acoustic  waves.  Both  bulk  acoustic  waves  and  surface  acoustic  waves 
(SAWs)  can  be  measured.  Brillouin  scattering  relies  on  thermally  activated  phonons,  meaning 
that  the  acoustic  amplitudes  are  small  and  that  measurements  are  time  consuming.  However  the 
exchanged  wavevector  is  fully  determined  by  the  scattering  geometry:  surface  phonons  can  be 
probed  at  wavelengths  down  to  a quarter  of  a micron.  Such  acoustic  wavelengths  are 
significantly  shorter  than  those  probed  by  any  other  measurement  technique.  Since  the  acoustic 
velocities  can  be  accurately  computed  as  functions  of  the  elastic  moduli  and  the  mass  density, 
their  measurement  provides  an  access  to  the  elastic  properties.  Brillouin  scattering  has  proved 
able  to  measure  the  elastic  properties  of  materials,  exploiting  either  bulk  acoustic  waves  [6]  or 
Surface  Brillouin  Scattering  (SBS),  i.e.  scattering  by  SAWs  [7-12],  In  particular  in  the  case  of 
layered  structures,  namely  of  a single  isotropic  layer  on  an  anisotropic  substrate,  the  elastic 
constants  of  films  have  been  measured,  for  film  thicknesses  ranging  from  hundreds  to  tens  of 
nanometers  [10-12],  The  measurement  of  film  properties  of  thickness  below  10  nm  remained  a 
challenge,  which  has  been  faced  in  this  work. 

EXPERIMENT 

The  tetrahedral  amorphous  carbon  (ta-C)  films  were  deposited  using  an  S-Bend  Filtered 
Cathodic  Vacuum  Arc  (FCVA)  with  an  integrated  off  plane  double  bend  (S-Bend)  magnetic 
filter.  The  deposition  chamber  was  evacuated  to  1 0‘8  torr  using  a turbomolecular  pump.  ta-C 
films  with  particle  area  coverage  of  less  than  0.01%  and  uniform  cross  section  are  consistently 
deposited  with  this  system  [13,14].  A series  of  films  was  grown  for  increasing  deposition  times 
defocusing  the  beam  in  order  to  achieve  a small  deposition  rate  and  so  to  produce  ultra-thin 
samples  in  a more  controlled  way.  Thickness  density  and  layering  of  these  ultra-thin  samples 
were  derived  by  X-Ray  reflectivity  [15],  Table  I.  Note  that  the  structure  of  the  ultra-thin  S-bend 
FCVA  films  resembles  that  of  thicker  S-bend  films  [14],  with  a scaling  of  the  bulk  layer 
thickness,  but  not  of  the  surface  and  interface  layers  (which  are  in  the  sub-nanometer  range  both 
in  thick  and  thin  films).  Four  samples  are  considered  here:  a bare  silicon  substrate,  with  a thin 
layer  of  native  oxide  due  to  exposition  to  air,  and  three  carbon  films. 

SAW  velocities  of  all  the  samples  were  measured  by  SBS:  the  specimen  is  illuminated  by 
laser  light  and  the  spectrum  of  the  scattered  light  is  analyzed.  Beside  the  strong  peak  of  the 
elastically  reflected  light  at  frequency  Q,  the  spectrum  contains  the  doublet  at  frequencies  Q+®, 
shifted  by  the  SAW  frequency  co.  The  SAW  wavevector  qy  is  selected  by  the  scattering  geometry 
(in  backscattering,  at  incidence  angle  9,  qy  = 2 (2tc/A)  sin  0,  where  X is  the  laser  wavelength)  and 
the  SAW  velocity  v = oV  qy  is  directly  obtained.  Measurements  were  performed  with  SAW 
propagation  along  the  [100]  direction  on  the  (001)  face  of  the  Si  substrate,  at  room  temperature, 
in  backscattering,  with  incidence  angle  9 from  30°  to  70°  (see  Fig.  1).  The  incident  light  is  p- 
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Table  I:  Samples’  properties 


Specimen 

Film  nature 

Film  density 
(g/cm3) 

Film  thickness 
(nm) 

A 

SiOx. 

— 

-2 

B 

ta-C 

2.8 

8 

C 

ta-C 

2.8 

4.5 

D 

ta-C 

3.2 

2.2 

polarised;  the  scattered  light  is  collected  without  polarisation  analysis,  and  analysed  by  a tandem 
3+3  pass  high  contrast  interferometer  of  the  Sandercock  type.  The  light  source  was  an  Argon  ion 
laser  operating  at  X = 514.5  nm.  The  incident  power  on  the  specimen  was  around  100  mW, 
focused  into  a spot  of  the  order  of  103  pm2;  irradiation  did  not  induce  film  modifications 

RESULTS 

Due  to  the  low  film  thickness  a single  SAW  is  present,  the  Rayleigh  wave  modified  by 
the  presence  of  the  film  (Modified  Rayleigh  Wave,  MRW).  Its  velocity  can  be  computed  [16]  as 
function  of  the  elastic  constants  and  mass  density  of  both  the  film  and  the  substrate,  of  film 
thickness  and  on  the  wavevector  qj.  It  is  thus  possible  to  solve  the  inverse  problem,  i.e.  the 
derivation  of  the  elastic  constants  from  measured  SAW  velocities  [17].  The  elastic  constants  and 
mass  density  of  the  anisotropic  Si  substrate  are  known  (Cu=166  GPa,  Ci2=63.9  GPa,  C<t4=79.6 
GPa,  p=2.33  g/cm3  [20]),  as  well  as  the  SAW  wavevector  and  the  film  thickness  and  density. 
Thus  the  acoustic  velocities  remain  functions  of  the  unknown  elastic  constants  of  the  film.  It  is 

140  . , 
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frequency  [GHz] 

Figure  1.  Portions  of  Brillouin  spectra  at  60°  incidence  showing  the  peak  due  to  modified 
Rayleigh  wave,  measured  under  identical  conditions  on  specimens  A (circles)  and  B (squares). 
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known  [11]  that  in  this  kind  of  films  the  properties  can  have  a gradient  in  proximity  of  the  outer 
surface  and  the  inner  interface  but,  to  avoid  exceeding  complexities,  the  films  are  modeled  as  a 
single  homogeneous  equivalent  film,  with  a sharp  interface  with  the  substrate. 

It  has  been  shown  [18]  that  among  the  various  elastic  moduli  that  can  represent  the 
stiffness  of  an  isotropic  material  the  (E,G)  pair  is  the  most  appropriate  for  the  solution  of  the 
inverse  problem.  The  velocities  vc(E,G)  are  computed  for  each  value  of  qy  (i.e.  for  each  incidence 
angle)  at  the  nodes  of  a rectangular  mesh.  The  couple  (E  ,G)fiim  is  obtained  by  a least  squares  fit 
of  the  computed  velocities  vc  to  the  measured  ones  vm  : 


<{E,G)-Vm 


= min, 


where  <Je'  are  the  variances  of  the  corresponding  vm‘  due  to  the  various  uncertainties  connected  to 
the  measurement  (geometrical  imperfections  [19],  statistical  errors).  The  (E,G)  couple  which 
minimizes  R is  the  most  probable  solution  of  the  inverse  problem,  and  the  regions  corresponding 
to  any  fixed  confidence  level  are  also  obtained. 

It  must  be  remembered  that  not  all  the  (E,G)  plane  is  physically  meaningful.  In  this  plane 
the  v=const  lines  are  straight  lines  through  the  origin,  and  the  v=0.5  and  v=0  line  delimit  the 
triangle  which  is  fully  meaningful.  The  v=0.5  line  is  the  limit  for  thermodynamic  stability.  When 
approaching  or  crossing  this  limit,  the  computed  SAW  velocities  remain  finite  and  continuous, 
but  the  bulk  modulus  B diverges.  This  means  that  if  the  confidence  region  falls  close  to  the 
stability  limit  even  a small  region  in  the  (E,G)  plane,  i.e.  a good  determination  of  both  E and  G, 
still  corresponds  to  a wide  interval  of  B,  which  remains  poorly  determined.  At  the  same  time  a 
physical  plausibility  limit  can  be  imposed,  considering  only  the  part  of  the  confidence  region 
which  corresponds  to  values  of  B below  a pre-determined  value.  In  this  work  the  diamond  value 
B = 445  GPa  is  taken  a physically  plausible  upper  limit.  On  the  other  side  the  v=0  line  does  not 
mark  a thermodynamic  instability  but,  since  negative  values  of  Poisson’s  ratio  can  at  most  be 
found  in  extremely  rare  and  peculiar  cases,  the  v=0  line  is  taken  as  a physically  plausible  limit. 
Figure  2 shows  the  confidence  regions  found  for  specimen  B.  The  confidence  regions  cross  the 
physically  acceptable  region,  meaning  that  v and  B remain  essentially  undetermined,  while 
reasonably  well  defined  intervals  are  found  for  E and  G:  E = 350  + 420  GPa  and  G = 130-5-210 
GPa.  For  specimen  C the  intervals  are  E = 220  + 260  GPa  and  G = 70  + 130  GPa,  v and  B 
remaining  also  undetermined.  This  is  in  agreement  with  the  outcome  of  a sensitivity  analysis 
study  [18]  which  showed  that  E and  G can  be  determined  much  better  than  v and  B.  The 
confidence  regions  must  be  seen  as  the  real  outcome  of  the  measurements.  Further  resolution 
within  the  confidence  regions  has  not  a real  significance.  The  obtained  values  of  the  elastic 
constants  allow  to  fit  the  measured  dispersion  relations,  as  shown  in  Fig.  3.  Figure  3 also 
explains  the  difficulty  in  obtaining  the  film  elastic  constants  for  specimens  A and  D,  whose 
measured  velocities  are  lower  than  those  of  a bare  silicon  substrate  of  nominal  properties  [20], 
The  fit  for  specimen  A is  in  fact  bad,  and  a better  fit  would  require  negligible  film  elastic 
constants  and  an  unreasonably  high  mass  density  of  silicon  oxide.  Also  specimen  D (carbon) 
seems  to  be  acoustically  slower,  i.e.  more  compliant,  than  silicon.  These  findings  suggest  that  in 
reality  the  actual  properties  of  the  silicon  substrate  are  slightly  different  from  the  accepted  values 
for  pure  silicon  [20].  The  stiffness  of  the  substrate  is  slightly  overestimated,  implying  that  the 
film  stiffness  is  underestimated.  Due  to  the  very  small  film  thickness  a very  small  overestimation 
of  substrate  stiffness  is  sufficient  to  lead  to  a significant  underestimation  of  film  stiffness. 
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Figure  2.  90%,  95%  and  99%  confidence  regions  in  the  (E,G)  plane  for  specimen  B 


DISCUSSION  AND  CONCLUSIONS 


It  has  been  shown  that  films  only  a very  few  nanometers  thick  exhibit  dispersion  relations  of  the 
modified  Rayleigh  wave  that  are  measurably  different.  The  analysis  of  these  dispersion  relations 
allows  to  measure  the  elastic  constants  of  ta-C,  films,  which  have  a significant  stiffness  also  in 
films  of  8 and  4 nanometres  thick.  For  thinner  films  the  analysis  shows  that  a very  precise 
characterization  of  the  substrate  is  needed.  The  results  for  the  thinnest  films  indicate  that  the 
values  of  the  elastic  constants  found  for  the  other  films  must  be  viewed  as  lower  bounds. 
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Figure  3.  Measured  velocities  for  the  four  specimens;  dispersion  relations  computed  with  the 
film  elastic  moduli  obtained  from  the  best  fit  procedure,  and  accepted  values  of  silicon  properties 
[20].  The  limit  at  null  wavevector  is  the  Rayleigh  velocity  of  a bare  silicon  substrate. 
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